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This paper presents a novel approach to wave propagation inside the Fabry-Perot framework. It 
states that the time-averaged Poynting vector modulus could be nonequivalent with the squared- 
field amplitude modulus. This fact permits the introduction of a new kind of nonlinear medium 
whose nonlinearity is proportional to the time-averaged Poynting vector modulus. Its transmittance 
is calculated and found to differ with that obtained for the Kerr medium, whose nonlinearity is 
proportional to the squared-field amplitude modulus. The latter emphasizes the nonequivalence of 
these magnitudes. A space-time symmetry analysis shows that the Poynting nonlinearity should be 
only possible in noncentrosymmetric materials. 
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I. INTRODUCTION 



A classical topic in electromagnetism is the study of 
wave transmission in a finite parallel-plane faces medium, 
know as Fabry-Perot resonator. When the medium 
presents nonlinear behavior, bistability appears Ij. To 
explain this phenomenon, the nonlinear Fabry-Perot res- 
onator (NLFP) was modelled by a third order suscep- 
tibility or Kerr-type nonlinearity 01 • At monochro- 
matic plane wave excitation, the NLFP stationary regime 
is summarized in a non-time-dependent nonlinear wave 
equation for the complex field amplitudes, the nonlinear 
Helmholtz equation (NLHE). Hence, the refiectance and 
transmittance problem reduces to finding the NLHE so- 
lution with appropriate boundary conditions for the field 
amplitude modulus and phase. 

The NLHE complexity led to approximate methods 
of resolution. Many approaches consider two counter- 
propagating waves in the medium and the analysis is 
done by separately considering the effects on each wave 

y, Ij, 1^ . Unfortunately, the linear superposition prin- 
ciple is no longer valid in nonlinear media and the sepa- 
ration of the electromagnetic field in these back and forth 
waves is meaningless. As a result, the NLHE separation 
into two equations, one for each wave, is only possible by 
neglecting various coupling nonlinear terms that would 
give an important contribution to the accuracy of the 
transmittance results. Moreover the Slowly Varying En- 
velope Approximation (SVEA) is often applied to these 
waves 0, y, Q when, yet within the counter-propagating 
wave approach, its validity was questioned p. Also, the 
boundary conditions were simplified rather than rigor- 
ously treating them 0, 0] ■ The above facts suggest that 
all these approximated approaches could not physically 
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be equivalents to the exact problem. 

The work done by Chen and Mills exactly solve the 
NLFP for an absorptionless Kerr-type medium 0. The 
proper of its resolution method was to assume a gen- 
eral complex field within the medium, disregarding the 
concept of counter-propagating waves. Chen and Mills 
derive a two coupled equation system for the field ampli- 
tude modulus and phase together with general boundary 
conditions, thus obtaining an analytic-transcendental so- 
lution for the transmittance of the NLFP. 

On the other hand, their work permitted us to note an 
implicit difference between the time-averaged Poynting 
vector modulus, i.e. the electromagnetic radiation inten- 
sity I, and the squared-field amplitude modulus (li^P) 
inside the nonlinear medium. If the nonequivalence of 
these magnitudes were true, it could change certain well- 
established fundamental concepts in classical electrody- 
namics. This fact motivated us to develop a novel ap- 
proach to wave propagation in nonlinear media inside 
the Fabry-Perot framework called S- Formalism. It in- 
troduces a new variable related to the time-averaged 
Poynting vector which states that its magnitude could 
be nonequivalent with the squared-field amplitude mod- 
ulus contrary to commonly accepted. Furthermore, the 
S-Formalism presents two important advantages: it per- 
mits to directly monitor the radiation intensity within 
the medium, and it avoids approximations, such as the 
SVEA, simplification of boundary conditions, and so on. 

The fact that the time-averaged Poynting vector mod- 
ulus be nonequivalent with the squared-field amplitude 
modulus, as the S-Formalism will show, implies that the 
nonlinearity of Kerr-type media is not proportional to the 
intensity which is contrary to what has been established 
to-date. This assertion leads to the following question re- 
garding the modelling of the NLFP: is it a Kerr-type non- 
linearity, or does it vary proportionately to the intensity? 
As this question does not have a definitive answer, the 
existence of the latter cannot be denied. Then, we define 
the Poynting medium as the medium where nonlinear- 
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FIG. 1: An harmonic plane wave strikes a nonlinear Fabry- 
Perot resonator, to be reflected and transmitted. The 
Regions I and III constitute, for simplicity, the same linear 
dielectric medium (e.g., air). 

ity is proportional to intensity. Thus, our objective is to 
solve the Poynting-NLFP through the S-Formalism com- 
paring the resultant transmittance with that obtained 
for a Kerr-NLFP to remark the noncquivalcnce between 
squared-field amplitude modulus and radiation intensity. 

In section II, we derive the S-Formalism in the follow- 
ing form: first the time-averaged Poynting vector assum- 
ing harmonic fields is calculated. Then, a dimension- 
less variable which is proportional to intensity is intro- 
duced. Thus, a set of field evolution differential equa- 
tions and the general boundary condition equations on 
the new field variable are obtained, which constitutes 
the S-Formalism. In section III, the new approach is ap- 
plied to derive the transmittance for the Poynting-NLFP 
and the results compared with that obtained for a Kerr- 
NLFP. A brief discussion about the possibily of existence 
and observation of Poynting media is given. Finally, in 
section IV, we conclude. 



II. THE S-FORMALISM APPROACH 

Refering to Figure ^ we start writing the linearly po- 
larized transversal harmonic electromagnetic fields of fre- 
quency oj as 

E,(z, <) = i (£;,-(z)e-'"* -t- c.c.) i, (la) 
H,(z, i) = i {Hnz)e-'^' + c.c.) j, (lb) 

where Ef{z) and Hf{z) are the non-timc-dcpcndcnt 
complex amplitudes for I =1, II, III. From region I, a 
plane wave of amplitude Eq and wave vector fco impinges 
perpendicularly on a nonmagnetic, isotropic and spatially 
nondispersive medium of thickness d (region II) . The op- 
tical field is assumed to maintain its polarization along 
this region so that a scalar approach is valid. The re- 
flected and transmitted plane waves have amplitudes tEq 
and IEq with r and t the complex reflection and trans- 
mission coefficients, respectively. Then, at region I and 



III the spatial-dependent complex amplitudes are given 

by 

E1^{z) = £:o (e'*^"^ -f re"*'^'"") , (2) 
E^jj{z) = i?o<e''=«^ (3) 

Similarly to Ref. at region II, we write down the 
following ansatz for the spatial-dependent complex am- 
plitude of the electric field: 



Et,{z) = E^£{z)e^t'^^\ 



(4) 



where the dimensionless amplitude S{z), and phase (/'(z), 
are both real functions of z. 

The time-averaged Poynting vector (E£(z, t) xH£(z, t)) 
can be easily calculated with the aid of Faraday's law, 
giving 



dz 



(5) 



where jUq is the vacuum permeability. From this expres- 
sion we calculate the intensity for the three regions: 
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{S)„ = lo fc„-if2(^)Mf) ^ i^s{z), (6b) 



dz 



{S) 
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(6c) 



where /q — k^E^ / {2^qIjS) is the incident intensity. In 
region II, Eq. (j6b|) defines the dimensionless field variable 



dz 



(7) 



directly related with the intensity inside the medium and 
which will characterize the S-Formalism. From Eq. {Tj) 
it is clear that if (j) is not a linear function on z, as is 
often happens in nonlinear media, then S and £^ are 
nonequivalents. 

The next step is to derive the NLHE in terms of the 
classical field variables (£, </)), and transform it into a set 
of equivalent equations in terms of (£, S). The NLHE is 
derived from the macroscopic Maxwell equations comple- 
mented by appropriate constitutive relations. We assume 
that the polarization P, and current density J, vary only 
in the electric field direction with frequency to, neglecting 
higher harmonics, and their spatial-dependent complex 
amplitudes satisfy the following constitutive relations: 

Pii{z) = eoXaen[z,Erj,Hyj]Erj{z), (8) 
J^j{z) = agen[z,Eyj,H'^'j]Eyj{z), (9) 

where eo is the vacuum permittivity and, Xgen and CTge„ 
are the generalized susceptibility and conductivity re- 
spectively, that are real and contain the linear as well as a 
possible nonlinear medium response. Note that the con- 
stitutive relations are not explicitly written since cases 
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could exist where it is not possible to describe the non- 
linear polarization and current density by the classical 
electric field power expansion. Thereby, the scalar NLHE 
is 



+ + Xgen) 



lUJflo a gen 



E^j{z)=0. (10) 



This equation constitutes the starting point to study sev- 
eral linear and nonlinear monochromatic wave propaga- 
tion phenomena within the Fabry-Perot framework. Sub- 
stituting Eq. Q into ifTUI) and by using Eq. jT)), we de- 
rive the following set of spatial evolution equations for 
the field variables £{z) and S{z): 



^ + ^0 ( il + Xaen[z,£,S])£- — ] =0,(lla) 



-T' + T-^^0 O-gen [z, £, S] f ^ = 0. 

dz kq 



£3 



(lib) 



To guarantee the physical content of the solution, these 
equations must be necessarily complemented with the fol- 
lowing boundary conditions: the continuity of the tan- 
gential components of the electric and magnetic field at 
the interfaces. The general boundary conditions were 
rigourously derived in Ref. four equations as func- 
tions of (f , 0) at z = and z — d which, by using Eq. 
(|7|) , are transformed into three equations in terms of the 
variables {£, S) to give 



the boundary conditions. Furthermore, through S{z) it 
is possible to monitor directly the intensity along the 
medium as a function of the spatial coordinate, as op- 
posed to using the conventional formalism. 

The S-Formalism is useful to analyze the linear case as 
well as the nonlinear one. Before to study the latter, i.e. 
the comparison between the Poynting and Kerr media in 
an effort to show the explicit difference between S (or /) 
and £^ (or |£'^|) in nonlinear media, we refer to the linear 
case. According to our analysis [8|, there are only two 
situations where the relationship / = cte\E'^\ holds true: 
firstly, a single plane wave propagates in a infinite or 
semi-infinite linear dielectric characterized by CTgen — 
and Xgen = X^'^'' where x^^'' is the linear susceptibility. 
Under these conditions, Eqs. (|ll|l relate the constants S 
and £ hy S = {1 + x^'^-^Y^^ £^; secondly, a single plane 
wave propagates in a semi-infinite linear absorber char- 
acterized by Xgen = X^^"^ a gen — <^ whcrc cr is the 
ohmic conductivity and such that S{z) oc £'^{z), being 
both proportional to a decreasing exponential function 
of z. On the contrary, when the medium is finite, e.g. a 
Fabry-Perot with boundary conditions at interfaces, S is 
no longer equivalent to £^ , not even for the linear dielec- 
tric case because 5" is a constant and 5^ is an oscillating 
function of z H|. 



III. THE POYNTING MEDIUM 



^d£ 

ko dz 



2=0 



S{d)-£^{d) = 0, 
d£ 

dz z=d 



0. 



4, (12a) 
(12b) 
(12c) 



From Eqs. © and (|12|l . the transmittance is obtained 
as 



\t?^S{d), 



(13) 



and the energy conservation is guaranteed thorough the 
expression: 



\r\' + W = l-{S{0)-S{d)). 



(14) 



This equation establishes that the refiectance and trans- 
mittance are limited by the boundary values of the time- 
averaged Poynting vector. For a nonabsorbent medium 
S{d) = 5(0), then |r|2 + |i|2 = l. 

Equations Hll|l and (|12|l represent the " S-Formalism" 
which were derived without assuming approximations 
such as the counter-propagating waves, SVEA, simpli- 
fications on the boundary conditions, and so on. Also, 
note that Eq. (|llb|) represents the time-averaged Poynt- 
ing Theorem applied to the problem of harmonic fields 
simplifying the interpretation of (Tgen as the dissipation 
properties of the medium. In particular, when Ggen = 0, 
the dimensionless intensity S* is a constant fixed by 



A. Constitutive relations and transmittance results 

At this point, we introduce the Poynting medium by 
the following constitutive relations: 



Xgen = X^^^ +7/o<S'(z), 
^ aen — Ot 



(15) 
(16) 



where 7 is the nonlinear coefficient. Eqs. (|llf) have a 
simple analytical solution given by 



S{z) = ^o, 




fc2 

1-t| cos [2ki{z-d)] + l + -^ 



(17) 
(18) 



where 7 > and k\ = /cq(1 + x(-^)+7/o<S'o). The constant 
Sq is fixed by 



1 - 



|)£^(0)+(2 



-i ) 5o - 4 = 0. (19) 



Combining Eqs. (|18|l and H19() . the transmitance can be 
expressed in a similar fashion as the linear Fabry-Perot 
resonator as 



T 



1 + Fsin^ {kid) 



(20) 
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X 



(1) 



1.25; and i = 2, ' = 5.25 



where F — kQ{l — fci/fco)^/(4fci). Carefully noting that 
Eq. (|20|l is a transcendental expression since ki depends 
on Sq. 

Now, we compare the transmittance results of the 
Poynting and Kerr media. The latter defined by 



Xgen = X^^' 



' gen 



= 0. 



(21) 
(22) 



The Kerr-NLFP transmittance results were taken from 
Ref. j|. Figure 2 shows T against the nonlinear param- 
eter 7/0 for two different values of x^^-*. Figures 2(0^) 
and 2(6i) correspond to the Poynting and Kerr medium, 
respectively. From these figures, it is apparent that the 
transmittance of the Poynting-NLFP as well as the Kerr- 
NLFP arc multistablc. However, for increasing values 
of x^^\ the peak transmittance separation diminishes 
for the Kerr medium while it increases for the Poynting 
medium. Also, the Kerr multistability appears for small- 
est values of the nonlinear parameter 7/9. The transmit- 
tance difference of both media emphasizes the / and 
noncqui valence. Figure 3 depicts T on the dimensionless 
thickness A:oc?(l -I- x'^'')^^^/(27r) enhancing the nonlinear- 
ity difference of Poynting and Kerr media. Note that the 
departure from Airy-type function for the Kerr medium 
is stronger than for the Poynting medium. 

On the other hand, Figure 4 shows that the Poynting 
nonlinear susceptibility, Xgen has a constant value 

along the medium (z coordinate). In return, the Kerr 
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FIG. 3: Transmittance against dimensionless thickness for 
(a) Poynting medium, (6) Kerr medium with x'^' = 5.25 and 
7J0 = 2. 
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FIG. 4: Nonlinear susceptibility against dimensionless spatial 
coordinate for each of the three solutions compatible with the 
boundary conditions. Continuous line: Poynting medium. 
Broken line: Kerr medium. The parameter values are 



FIG. 2: Transmittance against nonlinear parameter for (a^) 7/0 = 9, x 
Poynting medium, (bi) Kerr medium with ko d — 2n. For 



(1) 



1.25, and feo d, as defined in Fig. 2. 



nonlinear susceptibility varies periodically. This fact im- 
plies the formation of a phase grating in the Kerr medium 
on the contrary to the Poynting medium. Perhaps, this 
substantial difference can be measured, and this could be 
the starting point to experimentally identify a Poynting 
medium. 



B. Transformation properties under spatial 
inversion and time reversal 



It is a fact that several unusual types of nonlineari- 
ties were predicted before its experimental observation, 
as it was remarked, for example, in the pioneer works 
of Baranova et al. |^. With the aim to elucidate the 
isotropic medium requirements to observe the new phe- 
nomena, those authors pointed out the necessity of a 
analysis about transformation properties of electromag- 
netic quantities under rotations, spatial inversion and 
time reversal. Therefore, this symmetry analysis is also 
necessary to delimit the Poynting medium requirements. 

The magnitude that characterizes the electromagnetic 
response of a Poynting medium is their nonlinear suscep- 
tibility x*''^'' which is linear on the time-averaged Poynt- 
ing vector (S), as follows from the constitutive relations 
[Eq. IjlSf) ]. In a general form, it can be written as 



xr=7y-fe(S(r,t)),,, 
with i, j,k = x,y, z and 
1 /■*+^ 

(S(r,t)), = - [E(r,0 xH(r,t')],< 



(23) 



(24) 



where E and H are harmonics of period r — 2?:/ to and 
time interval T ^ t. The susceptibility tensor trans- 
forms as even under spacial inversion (r — > — r) and time 
reversal (t — > —t), contrary to the Poynting vector and 
its time-averaged which transforms odd under spatial in- 
version and time reversal, i.e. (S(r,t)) — *■ — (S(— r, f)) 
and (S(r,t)) -(S(r, -t)), respectively jlfl. Then, a 
medium possessing a linear connection between Xij smd 
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{S)z should be noninvariant with respect to spatial inver- 
sion and time reversal. Otherwise, the space-time sym- 
metry will be violated in the constitutive relation [Eq. 
(23)]. 

The lack of parity symmetry under inversion of coor- 
dinates is proper of materials without inversion center, 
i. e. Poynting nonlinearity should be only possible in 
noncentro symmetric materials. There are several mate- 
rials candidates to possess a Poynting nonlinearity, as 
for example the cubic crystals with zincblended struc- 
ture such as GaAs, InSb and others. In these materials 
intensity-dependent transmission and bistability was ex- 
perimentally observed (llj. Also, isotropic homogeneous 
liquids formed by nonracemic mixtures or solutions of 
mirror-asymmetric (chiral) molecules with strong nonlin- 
ear optical susceptibility, as product of several nonlinear 
processes , are also feasible of posses a Poynting non- 
linearity. In addition, parity under time reversal should 
be violated in Poynting media. This means that any weak 
dissipative process, that converts field energy into heat, is 
necessary to remove the rule relating to the t —t trans- 
formation. For example, either very weak absorbtion or 
current flow by external quasi-static field, that basically 
do not affect the wave propagation at light frequency to, 
would ensure the medium non-invariance under time re- 
versal. 

We believe that, in spite of experimental works are 
required, the above preliminary analysis could stimulate 
further discussions regarding the existence of Poynting 
media. 



propagating in a linear infinite or semi-infinite medium. 
Otherwise, they are nonequivalents. Besides, the S- 
Formalism presents two important advantages: it per- 
mits to directly monitor the time-averaged Poynting vec- 
tor in the medium and it avoids approximations, such 
as SVEA, simplification of the boundary conditions, and 
so on. To emphasize this nonequivalence we introduce 
the Poynting medium, whose nonlinearity is proportional 
to the intensity instead of the electric squared-field am- 
plitude modulus such as in the Kerr medium. We find 
marked disagreement in the transmittance of both media, 
which support the differences between I and Also, 
a space-time symmetry analysis shows that the Poynting 
nonlinearity should be only possible in noncentrosym- 
metric materials. 

The statements and analysis pointed out here consti- 
tute an advance on theoretical views of basic concepts 
in electrodynamics. The S-Formalism could be impor- 
tant in problems where the time-averaged Poynting vec- 
tor must be rigourously monitored like in photoconductor 
or photorefractive materials. Further to this particular 
case studied here, this new approach leaves open the pos- 
sibility of new physical results in actual topics on non- 
linear wave propagation such as spatial solitons, wave 
mixing and others. Finally, we leave open the possibility 
that experimental techniques, based on intensity depen- 
dent phase changes of a Gaussian beam such as Z-Scan 
Technique 13], could not truly measure Kerr-type non- 
linearity. On the contrary, they could be measuring a 
Poynting-type instead. 



IV. CONCLUSIONS 

In summary, we derive a new formalism in terms of 
dimensionless variables related with the time-averaged 
Poynting vector and field amplitude modulus within the 
Fabry-Pcrot framework. The S-Formalism shows explic- 
itly that the energy intensity and squared-field ampli- 
tude modulus are only equivalents for a single plane wave 
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